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Recently, neutron scattering spin echo measurements have provided high resolution data on the 
temperature dependence ol the linewidth T(q, T) of acoustic phonons in conventional superconduc- 
tors Pb and Nb. 1 At low temperatures the merging of the 2A(T) structure in the linewidth with a 
peak associated with a low lying h\Uq KA Kohn anomaly suggested a coincidence between 2A(0) and 
fia)q KA in Pb and Nb. Here we carry out a standard BCS calculation of the phonon linewidth to 
examine its temperature evolution and explore how close 2A(0)/fio;q KA must be to unity in order 
to be consistent with the neutron data. 



PACS numbers: 

I. INTRODUCTION 

Using resonant spin echo neutron scattering tech- 
niques, Aynajian et at 1 have recently measured the 
linewidth of transverse acoustic phonons in high purity 
single crystals of Pb and Nb. At low temperatures, which 
are however above the superconducting transition tem- 
perature T c , a plot of the phonon linewidth T(q, T) as a 
function of the phonon wavevector q exhibits peaks which 
arise from Kohn anomalies. 2 When the temperature de- 
creases below T c and the superconducting gap opens, one 
sees an expected decrease in the linewidth T(q, T) for 
phonons having energy fiw q less than twice the super- 
conducting gap A(T). As ficj q approaches 2A(T), there 
is a rapid increase in T(q, T) associated with the peak in 
the quasiparticle density of states at the gap edge and 
the fact that the BCS coherence factor for a phonon to 
break a Cooper pair and decay into two quasiparticlcs 
approaches 1 at threshold. 3 However, as Aynajian et al. 
note, what is surprising is that as T goes to zero, the 
feature in r(q, T) that is associated with fiw q = 2A(T) 
appears to merge with a Kohn anomaly peak. This be- 
havior is seen in both Pb and Nb, posing the question of 
why should the energy of a transverse acoustic phonon 
associated with a normal state Kohn anomaly coincide 
with twice the limit of the low temperature supercon- 
ducting gap 2A(0)? 

Motivated by this experimental result, we have carried 
out a standard BCS calculation of the temperature de- 
pendence of the transverse acoustic line width and exam- 
ined what happens if 2A(0) is near the energy associated 
with a normal state Kohn anomaly in T(q, T). In partic- 
ular, we are interested in the evolution of T(q, T) as the 
temperature is lowered and 2A(T) approaches the energy 
of the Kohn anomaly huj <lKA . How close to hujq KA docs 
the low temperature limit of 2A(T) need to be for it to 
appear that the 2A(0) structure in T(q, T) merges with 
the Kohn anomaly structure as T goes to zero? 



II. FORMALISM 

We begin by first examining the matrix elements for 
the electron coupling to the transverse acoustic modes. 
In clean materials, the coupling of the electrons to the low 
frequency transverse phonons occurs through Umklapp 
scattering processes. As one knows, this is because the 
polarization e\ (q) of a transverse phonon is orthogonal to 
q. In Fig. la,b we show Fermi surface sections for Pb and 
Nb, respectively, obtained from Density Functional The- 
ory (DFT) calculations (ABINIT). 4 In the top figure for 
Pb, an Umklapp scattering process is shown in which an 
electron is scattered from k to k' = k + K„ + q with q the 
wavevector of the transverse phonon and K„ a reciprocal 
lattice vector. In this case, the phonon wavevector c[ka 
that is shown connects two parts of the Fermi surface 
that have parallel tangents leading to a Kohn anomaly 
in the scattering rate and the phonon linewidth. A sim- 
ilar process for Nb is illustrated in the lower part of Fig. 
1. 

In the following calculations, we use an electron- 
phonon vertex g\(k, k') 

o A (k, k') = -u / % ^e A (k - k') • 

tfM ' ' y 2MNuj\(k — k ) M ; 

]T (k + K m - k' - K n )al m (k)a Kn (k') x 

K m ,K„ 

(k + K m |[/|k' + K n ) (1) 

for a transverse acoustic mode A which has a frequency 
hu>\(q) and a polarization vector e\(q). Here M is the ion 
mass, U is the lattice pseudo-potential, and N is the num- 
ber of lattice sites. As discussed, for transverse phonons 
one needs an Umklapp process to couple the electrons to 
the direction of the ionic vibration given by c\—t such 
that qt oc (K + q) • et = K • £t- Then, as we will see, the 
linewidth of the transverse phonons will exhibit a peak 
as q approaches the Kohn anomaly wavevector qKA- 
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FIG. 1: (Color online) The Fermi surfaces of (a) Pb and (b) 
Nb. The short black arrows indicate the qxA wavevectors 
(r.l.u.) q KA = 0.295 in Pb and q KA = 0.196 in Nb. The long 
arrow indicates the reciprocal lattice vector K n associated 
with the Umklapp scattering process. The location of the 
Bragg planes are indicated by thin dashed lines. 



To capture the essence of the Kohn-Umklapp scatter- 
ing, we first consider the expression for the transverse 
acoustic phonon linewidth T(q, T) in the normal state 
for the case in which the Fermi surface spanned by q KA 
is approximated by a cylinder of radius kp (Fig. 2a). In 
this case, T(q, T) is given by 



T(q,T) 



ngKi 

N 



k 

xS(ojq - e k+q + e k ) 



(2) 



with / the Fermi factor and e k the electronic band dis- 
persion. From here on we choose % = 1. For simplic- 
ity, we have set g,\=T(k,k') = ^k, the phonon mode 
energy cj q = Cx|q|, with ct the transverse speed of 
sound, and assumed a simple 2D free-electron dispersion 
e k = k 2 /2m — /i. Taking the T = limit and making the 



10 

8 
6 
4 
2 



CT 

Z< 3 



.(a) j 










(b) , 











0.25 



0.5 



0.75 

q/q 



1.25 



KA 



FIG. 2: The zero temperature phonon linewidth of the trans- 
verse acoustic branch in the normal state with ct/vf = 0.01 
evaluated for free electrons with (a) a cylindrical Fermi sur- 
face and (b) concave Fermi surface. 



change of variables x = k/kp, Eq. (2) reduces to 

r(q,r = o) = mfcF|gK|2 fxdx jd<p (3) 

4tt<7 J J 
o o 

x [S (a_(q) — xcos((/>)) — S (a+(q) — xcos(4>))] 

where a±(q) = ^ ± and kp and Vp are the Fermi 
momentum and velocity, respectively. After a little alge- 
bra we then obtain 



r(q,T = 0) = N F \g K \ 



1 2fc; 



l-cd(q)e(l-a 2 _(q)) 



-Wl-^(q)e(l-<(q)) 



(4) 



where Np is the single-particle density of states per spin 
at the Fermi level and Q(x) is the usual step function. 

T(q,T = 0) is plotted in Fig. 2a for c T /v F = 0.01. 
While the overall magnitude of the linewidth is deter- 
mined by the ratio of ct/vf, the momentum dependence 
comes from simple phase space considerations. One can 
see that the phonon linewidth grows rapidly for mo- 
mentum transfers approaching 2k p in this example and 
quickly falls to zero for larger momentum transfers as no 
phase space is available for scattering. This strong en- 
hancement of the phonon linewidth in the normal state 
at q corresponding to the Kohn anomaly will also be 
present in the superconducting state, with an additional 
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FIG. 3: (Color online) A schematic of the kinematic con- 
straint for the decay of an acoustic phonon in the supercon- 
ducting state. 



kinematic constraint imposed by the breaking of Cooper 
pairs. 



The kinematic constraint of phonon decay in the super- 
conducting state brings the energy scale 2A directly into 
play. This is shown in Fig. 3, which sketches quasiparti- 
cle scattering across the gap edge 2 A. Due to the disper- 
sion of the phonon, vertical scattering processes having 
no net wavevector transfers are kinematically forbidden. 
In order to bridge the gap the energy of the phonon must 
be at least 2 A. In other words, a finite wavevector trans- 
fer must occur where q = 2A/cy. In addition, the domi- 
nant Kohn-Umklapp process k' — k = K„ + q KA involves 
a momentum transfer of q^^, which in this sketch is 2kp. 
Thus we have two conditions that lead to the conclusion 
that when w(q) of the transverse acoustic phonon branch 
equals twice the superconducting gap, or in other words 
when q = qKA = 2A/c-r, an enhancement of the phonon 
decay will occur. 



We next consider the transverse acoustic phonon decay rate in the superconducting state for this cylindrical model 
using conventional BSC theory. 5 The phonon self-energy can be obtained by evaluating the electron-hole bubble. In 
the superconducting state the phonon self-energy II (q, iu) m ) is then given by 



n(q,iu m ) = — Tr^ \g K \ 2 T 3 G(k,iui n )G(k + q,iuj n + iu m )f 3 

" n,k 



(5) 



where u) n = (2n + l)7r//3 and uj m — 2nnr//3 are Fermion and Boson Matsubara frequencies, Tr denotes the trace, and 
G is the electron propagator 



G(k,iu n ) = 



iui n T + e k f 3 + A k fi 



(6) 



Here t% are the usual Pauli matrices and E^ = yi^-j- A£ is the quasiparticle energy. After analytic continuation, the 
phonon self-energy is given by 



n(q,^ q ) - -^^| ffK | 2 {^ + (k,q)[/(i? k )-/(£; k 



2N 

+ A_(k,q)[/(- J E k )-/(£ k+q )] 
with the coherence factors defined as 



flWq — E k + -E k+q + iS tujq + E k — -E k+q + iS 



k+q 



1 

iS fiu n — E\r — E] 



k+q 



id 



A±(k,q) = l± 



£k£k+q — A k A k+q 

EkEk + „ 



(7) 



(8) 



The q-dependent phonon lincwidth T(q, T) is then determi ned from the imaginary part of n(q, ui q ) 



The first two terms in Eq. 7 describe quasiparticle scat- 
tering processes. For these processes, the BCS coherence 
factor A + vanishes at the threshold and this, along with 
the depletion of the thermal quasiparticle populations 
as the gap opens, suppresses their contribution to the 
phonon linewidth. The fourth term in Eq. 7 corresponds 
to a process in which a phonon breaks a pair, creating 



two quasiparticles with wavevectors k + q and — k. This 
requires that the phonon energy w q be greater or equal 
to 2A(T). In this case, the BCS coherence factor A_ 
goes to 1 at threshold where E(k + q) = E(k) = A(T) 
and there is a sudden increase in the linewidth. 

Before turning to the results for the linewidth in Pb 
we first consider two simplified cases at T = 0, shown 
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in Fig. 4. Here we have set the phonon energy w q = 
crq and plotted T(q, T = 0) versus q/qKA- The A = 
curve is identical to Fig. 2. As the superconducting gap 
opens r(q, T — 0) is suppressed for crq < 2A(0) due 
to the loss of phase space for electron-phonon scattering. 
This produces an onset (or "knee") in r(q, T — 0) at 
an energy corresponding to the gap edge. Note that in 
this case one expects a knee rather than a peak because 
qt; ~ vf/ct l. 7 (The knee is also somewhat smeared 
here due to the finite broadening S — ct /40 used.) 

The height of the onset is controlled by the momen- 
tum q for breaking a Cooper pair into two quasiparticles 
carrying momenta k and k — q respectively. 6 ' 7 As A is 
made larger, the onset at orq — 2A associated with pair- 
breaking in the superconducting state moves out towards 
the Kohn anomaly at qxA- For a cylindrical Fermi sur- 
face, the Kohn anomaly occurs at qxA = Zkp and when 
2A(0) = CTqKA-, the pair-breaking onset coincides with 
the Kohn anomaly peak. If 2A(0) exceeds cxqKA, the 
Kohn anomaly peak is suppressed by kinematics as the 
energy to break a pair is greater than crqKA- 

For a cylindrical Fermi surface with 2A(0) > orqKA, 
T(q) is suppressed due to phase space considerations pre- 
viously discussed for the normal state. However, if the 
Fermi surface has some degree of curvature along the k z 
direction such a sharp cut-off will not occur. To illustrate 
this, in Fig. 4b we plot T(q, T = 0) for a Fermi surface 
that has a concave warping along the z direction (see Fig. 
2b). The electronic band dispersion has again been mod- 
eled by a free electron dispersion but with m x = m y = m 
and m z — —5m. 8 For such a dispersion qxA corre- 
sponds to the spanning condition across the narrowest 
portion of the Fermi surface (k z — 0). As can be seen in 
Fig. 4b, the concave curvature of the Fermi surface pro- 
vides phase space for scattering with momentum trans- 
fers q > q_ K A an d the sharp cutoff in r(q, T — 0) is 
no longer present. With the opening of the supercon- 
ducting gap, T(q, T = 0) is suppressed for crq < 2A(0), 
just as in the previous case. For 2A(0) = crqKA a rem- 
nant of the Umklapp-Kohn peak remains. As the gap is 
increased further (2A(0) > CtQka) the phase space as- 
sociated with the Kohn peak is gapped out and the peak 
in r(q, T = 0) is thus suppressed. 



III. RESULTS FOR LEAD 

With the simple examples of the previous sections we 
are now ready to turn to the phonon linewidth in Pb. To 
obtain the electron dispersion the DFT bandstructure for 
Pb was calculated on a regular grid of 100 x 100 x 100 
momentum points per quadrant of the first Brillouin zone 
and a linear intcrpolant was used to obtain energies at 
intermediate momenta. For the phonon dispersion we 
again assume a linear phonon dispersion w(q) = Or\q\, 
with ct — 7.9 meV/[r.l.u.]. The transition temperature 
T c = 7.2 K sets the temperature scale and we use an 
intrinsic broadening 8 = 0.01 meV throughout. Finally, 




FIG. 4: (Color online) The normalized transverse acoustic 
phonon linewidth in the superconducting state at T — for 
various values of the superconducting gap 2A = Ac T qK a- (a) 
r(<?) for a perfectly cylindrical Fermi surface, (b) T(q) for a 
cylindrical Fermi surface with a concave warping along the k z 
direction (see text). 



we note that an explicit evaluation of the matrix element 
for Umklapp scattering <7k,k' given by Eq. (1) adds an 
additional level of difficulty to the problem. Therefore, 
for simplicity, we approximate the matrix element with 
a constant c/k and restrict the momentum sum to the 
region near the orange (light) Fermi surfaces shown in 
Fig. la. 

The results are shown in Fig. 5 as a function of temper- 
ature for gap sizes ranging from 2A(T = 0) = 0.8cxqKA 
to l.lcTqKA- The qualitative behavior of T(q, T) is sim- 
ilar to that which was found for the simplified models 
considered in the previous section. Above T c the phonon 
linewidth is finite for all values of q = (q, q, 0) and has a 
peak at q = qxA = 0.285 [r.l.u.], which is associated with 
the Kohn anomaly indicated in Fig. 1. As the tempera- 
ture is lowered across T c , the gap opens following an as- 
sumed BCS temperature dependence. For cxq < 2A(T), 
T(q, T) is suppressed and the expected 2 A onset (knee) 
forms. (Here, T(q) has a finite value for crq < 2A(T) 
which is exponentially suppressed as T is lowered. This 
is due to the non-zero contributions of the first two terms 
in Eq. (7) and corresponds to the thermal occupation of 
quasiparticle states across the gap edge.) As T is low- 
ered further 2A(T) grows and the knee in T(q, T) moves 
towards the Kohn peak. If 2A(0) is smaller than crqKA 
this knee stops short of the peak at the lowest tempera- 
tures (Figs. 5a) while for for 2A(0) = crqKA it merges 
with the peak (Fig. 5b). Finally, if 2A(0) > c T q KA (Fig. 
5c), then for sufficiently low temperatures the Kohn peak 
is suppressed similar to the results shown in Fig. 4b. 

Thus within a BCS framework, T(q, T) depends upon 
the shape of the Fermi surface and Qka, the velocity 
of sound for the transverse acoustic branch, and the 
magnitude of the superconducting gap. The appropri- 
ate parameters for Nb and Pb are summarized in Tbl. 
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FIG. 5: (Color online) The linewidth of a transverse acoustic 
phonon in Pb as a function of reduced temperature t = T/T c 
for various values of 2A(T = 0) as indicated. 



I. For Pb we estimate 2A(0) = 0.95ctQka, which cor- 
responds closest to Fig. 5b, while for Nb we estimate 
2A(0) = 0.8cxqKA corresponding to Fig. 5a. 

Comparing our results to Figs. 3 and 4 of Ref. 1 we 
find that agreement with the experimental data for Pb is 
good while the agreement for the case of Nb is less clear. 
For Pb we find 2A(0) ~ 0.95ctQka and we therefore ex- 
pect a knee to form in T(q, T) which tracks out to the 
Kohn peak as the temperature is lowered. This behav- 
ior is similar to what is observed experimentally (Fig. 3 
of Ref. 1). In the case of Nb 2A(0) - 0.80c T q KA and 
we therefore expect the knee to approach the Kohn peak 
but stopping short at the lowest temperatures leaving a 
pronounced knee in the observed linewidth. Examining 
Fig. 4b of Ref. 1 it is difficult to determine if such a 
knee is present in the data. Finally, we note that our 
calculations predict that the Kohn peak should be sup- 
pressed when 2A(0) > ctQka- Therefore, one clear way 
to test the conclusions of this work would be to exam- 
ine the linewidth of the transverse acoustic branch in a 



TABLE I: The relevant parameters for the elemental super- 
conductors Pb and Nb. The values for Pb have been esti- 
mated from Ref. 1. The gap for Nb was obtained from Ref. 
11. The transverse speed of sound in Nb was obtained from 
Ref. 10. 



material where CtQka < 2A(0). 



IV. CONCLUSIONS 

We have seen that the momentum and temperature 
dependence of the transverse acoustic phonon linewidth 
T(q, T) in the superconducting state depends on Wq KA 
and 2A(T). While both of these energies depend upon 
the bandstructure and phonon dispersion, there is noth- 
ing that should lock them together in the traditional the- 
ory Thus while it is known that the Kohn anomaly 
wavevector (\ka shown in Fig. 1 gives rise to a small 
kink in a 2 F{oS) associated with the energy ojq KA at which 
the transverse phonon begins to contribute to the pair- 
ing interaction, this is a small feature and plays no role 
in determining A(0). 5 Therefore, within the BCS frame- 
work, the fact that 2A(0) is close to the energy of a Kohn 
anomaly w qKA , must be viewed as a coincidence. Fur- 
thermore, as noted, the fact that the wavelength of the 
phonon is small compared to the coherence length leads 
to a knee-like feature at 2A(0) rather than a peak. There- 
fore if 0Jq KA > 2A(0), the Kohn anomaly remains as 
the dominant feature at low temperatures. However, as 
shown in Fig. 4b and 5a, if 2A(0) is slightly less than the 
Kohn anomaly phonon energy ojq KA , the 2A(T) struc- 
ture can appear to merge with the Kohn anomaly peak 
in r(q, T) as T goes to zero. Thus we would conclude 
that it is an interesting coincidence that 2A(0) is only 
slightly smaller than the energies of the Kohn anoma- 
lies in both Pb and Nb, but it does not necessarily mean 
that the superconducting gap is determined by the Kohn 
anomaly itself and does not force 2A(0) = w qKA . 
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